Functional modularity of background activities in normal and epileptic brain networks 
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We analyze the connectivity structure of weighted brain networks extracted from spontaneous 
magnetoencephalographic (MEG) signals of healthy subjects and epileptic patients (suffering from 
absence seizures) recorded at rest. We find that, for the activities in the 5-14 Hz range, healthy brains 
exhibit a sparse connectivity, whereas the brain networks of patients display a rich connectivity 
with clear modular structure. Our results suggest that modularity plays a key role in the functional 
organization of brain areas during normal and pathological neural activities at rest. 
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From the brain to the Internet and to social groups, 
the characterization of the connectivity patterns of com- 
plex systems has revealed a wiring organization that can 
be captured neither by regular lattices, nor by random 
graphs [l|. In neurosciences, it is widely acknowledged 
that the emergence of several pathological states is ac- 
companied by changes in brain connectivity patterns [2^] . 
Recently, it has been found that functional connectivity 
patterns obtained from magnetoencephalography (MEG) 
and electroencephalography (EEG) signals during dif- 
ferent pathological and cognitive brain states (including 
epilepsy) display small- world (SW) properties [3]. Em- 
pirical studies have also lead to the hypothesis that brain 
functions rely on the coordination of a scattered mosaic 
of functionally specialized brain regions (modules), form- 
ing a web- like structure of neural assemblies [2] . Modu- 
larity is a key concept in complex networks from RNA 
structures to social networks ^, |5|] . A module is usually 
defined as a subset of units within a network, such that 
connections between them are denser than connections 
with the rest of the network. In biological systems, it is 
generally acknowledged that modularity results from evo- 
lutionary constraints and plays a key role in robustness, 
flexibility and stability [6]. 

Absence seizures are the most characteristic expres- 
sion of non- convulsive generalized epilepsy. Their main 
characteristic in brain signals is the occurrence of high 
amplitude, rhythmic spike wave discharges synchronized 
over wide cortical areas, which manifests suddenly from 
a normal background. Current studies of brain connec- 
tivity mainly focus on the onset and evolution of epilep- 
tic discharges. Nevertheless, little is known about the 
organization of brain networks during spontaneous back- 
ground activities of patients, and the role of these connec- 
tivity patterns on the emergence of absence seizures. In 
this Letter we study the modular organization of brain 
networks extracted from spontaneous MEG signals of 
epileptic patients and healthy subjects. The results of 
our analysis reveal a non-random structural organization 
in both normal and pathological brain networks. In par- 



ticular, the functional networks of control subjects are 
characterized by a sparse connectivity between the mod- 
ules. In contrast, brain connectivity of epileptic patients 
(recorded out of seizures) displays a configuration where 
nodes in a functional module are connected to different 
functional modules. This modular configuration might 
play a key role in the integration of large scale brain ac- 
tivities, facilitating the emergence of epileptic discharges. 

The data used in this study were acquired from 5 
healthy subjects and 5 epileptic patients suffering from 
absence seizures. The study was performed with writ- 
ten consent of the subjects and with the approval of 
the local ethics committee. During the recordings, sub- 
jects and patients were instructed to rest quietly, but 
alert, and keep their eyes closed. The brain signals 
were acquired with a whole- head MEG system (151 sen- 
sors; VSM MedTech, Coquitlam, BC, Canada), digitized 
at 1.25 kHz with a bandpass of - 200 Hz. All the 
analyses were performed on 338 non-overlapping quasi- 
stationary segments (206 for all the patients and 132 for 
the healthy group) of 5 seconds without eyes or muscu- 
lar artefacts, nor epileptic activities (as, e.g., seizures or 
epileptic-like activity) and far (at least 10 s) from recent 
epileptic discharge. In agreement with previous findings, 
surrogate data tests revealed that less than 4 % of in- 
terdependencies between the spontaneous brain activi- 
ties were nonlinear [7]. Thus, weighted brain networks 
were constructed by means of a definition of functional 
links based on linear coherence. The squared modulus 
of the coherence between two time series Xi {t) and Xj (t) 
(normalized to zero mean and unit variance) was defined 

as: |r^j(/)/ = -g^^^jyp^-^^ where Sa and Sij (the spec- 
tral and cross-spectral densities) were estimated using 
the Welch's averaged periodogram method |8|]. 

Recent results show that correlations between mag- 
netic fields sensors located at a distance less than 4 cm 
can not distinguish between spontaneous activities of 
epileptic patients and control subjects [9^. To reduce the 
influence of these spurious correlations between MEG sig- 
nals, we have excluded the nearest sensors (separated less 
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than 5 cm) from the computation of coherence values. To 
perform the statistical analysis of coherence values, we 
used Fisher's Z transform of Tij: Zij = 0.5 In (^j^^^ 
Under the hypothesis of independence, Zij has a normal 
distribution with expected value and variance l/2A/'5, 
where TV^ is the number of non-overlapping blocks used 
in the estimation of spectral densities [8]. To correct for 
multiple testing, the False Discovery Rate (FDR) method 
was applied to each matrix of Tij values [lOj. With this 
approach, the threshold of significance Tth was set such 
that the expected fraction of false positives is restricted 
to q < 0.01. Finally, in the network construction, a func- 
tional connection between two of the N = 151 nodes 
(brain sites) was assumed as an undirected weighted link, 
i.e. we set the weight of the link between i and j as 

and Wii =0 otherwise. 



13 = ^ji = Tij when F^^- > Tth, 
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FIG. 1. Topographic distribution of node strength Si (a), 
weighted clustering coefficient Ci (b), and efficiency Ei (c). 
Control subjects (CTL), patients (PAT), and difference maps 
(p- values) of their comparison are reported for real and equiv- 
alent random configurations. 

To characterize the network structure of healthy sub- 
jects and epilectic patients, we evaluated a list of mea- 
sures for weighted graphs [1| . Namely, the node strength 
Si = J2k^ik, the weighted clustering coefficient q [1[, 
the efficiency Ei of a node i [ijj; as well as the re- 
spective averages over the graph nodes S = ^/N^-Si, 
Cw = 1/N Ci and the global efficiency of the graph E. 
As null model, brain networks were compared to equiva- 
lent random graphs, obtained by randomly rewiring the 
links of the original networks while keeping the same de- 
gree distribution [l2|]. The statistical significance of a 
given metric is assessed by quantifying its statistical 
deviation from values obtained in the ensemble of ran- 
domized networks. Let ja and a be the mean and SD of 
the parameter computed from such an ensemble. The 
significance is given by the ratio T^o = \0 — fi\/a whose 



p- value is given by the Chebyshev's inequality [iS^ . To as- 
sess significant differences between topological features of 
the two groups, we used standard nonpar ametric permu- 
tation methods, which account for spatial dependences 
in the data [1^. We used exhaustive permutations (10^) 
to estimate the empirical distribution under the null hy- 
pothesis of no difference between the two groups. 



Cuj (mc) 



Q {^q) Nrr 



CTL 4.04 0.233 (0.096) 0.649 (0.892) 0.538 (0.001) 13.7 
PAT 7.34 0.300 (0.160) 0.893 (1.101) 0.503 (0.006) 8.30 

TABLE L Network properties in control subjects (CTL) and 
patients (PAT): mean node strength 5*, average weighted clus- 
tering coefficient Cw^ global efficiency E^ maximal modularity 
Q and number of modules Nm • Me denotes the average of met- 
ric obtained from 20 random graphs. 

Although we applied our approach to connectivity 
graphs obtained from brain oscillations at the conven- 
tional frequency bands (/ < 5 Hz, 5 < / < 15 Hz, 
15 < / < 24 Hz, 24 < / < 35 and / > 35 Hz), statisti- 
cally significant differences between normal and epileptic 
brain networks were observed only for the brain activi- 
ties in the extended alpha range (5 — 14 Hz). Henceforth, 
all results presented here refer to functional networks ob- 
tained at this frequency band. The basic network proper- 
ties (averaged over subjects) are summarized in Table HI 
The structure of functional brain networks was found 
to be significatively different from that of randomized 
counterparts. Namely, brain networks of both patients 
and healthy subjects yielded a clustering coefficient Cw 
larger that that of randomized graphs {p < 10~^), and 
values of efficiency as large as those of random graphs, 
indicating a small- world behavior. These results agree 
with previous findings suggesting that brain sites have 
an optimal interaction with most other brain regions [ll . 
Furthermore, average network properties (5*, and E) 
of epileptic patients were found to be different (p < 10~^) 
from those of control subjects, supporting the hypothesis 
that neural disturbances are correlated with changes in 
functional network architectural features [2]. 

A more detailed information on the differences between 
patients and control subjects can be acquired by the anal- 
ysis of the network at the level of node properties. In 
Fig. [1] we report the spatial distribution of node mea- 
sures Ci and Ei (averaged over all control subjects and 
over all patients) for each sensor of the network. Results 
indicate that epilectic patients have a richer node con- 
nectivity than control subjects. Difference maps clearly 
identify the centro-parietal regions as those brain areas 
with the highest contrast between patients and control 
group. Although the node strength of epileptic patients 
is twice as large as that of healthy subjects, if links of real 
networks are randomly rewired by keeping the same de- 
gree distribution, no significant differences are observed 
between the two groups. This rules out the possibility 
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that the difference in the number of connections alone 
could account for differences in and E. No signif- 
icant difference between healthy subjects and patients 
was observed by direct comparison of the spatial maps of 
the a-activity power [22 ]. 

A potential modularity of brain networks is suggested 
by the fact that the networks display a clustering coef- 
ficient larger than that obtained in random graphs [15]. 
Previous studies over brain networks have used cluster- 
ing methods to identify similar groups of brain activities. 
However, classical approaches such as those based on 
principal components analysis (PC A) and independent 
components analysis (ICA), make very strong statistical 
assumptions (orthogonality and statistical independence 
of the retrieved components, respectively) with no phys- 
iological justification [16]. To find the network modules 
we have instead used an algorithm based on a spectral 
embedding of graphs [17] (similar theoretical frameworks 
have been recently proposed [HI). The algorithm is 
based on the definition of a Markov chain with a tran- 
sition probability matrix P with entries pij = If 
is the t^^ iterate of matrix P, the element Pij{t) en- 
codes the probability of moving from node i to node j 
through a random walk of length t. For an undirected and 
connected graph, the Perron Frobenius theorem assures 
that /i* = is the unique stationary distribution of 

the Markov chain, such that lim^^oo ^iPij{^)l^i{^) = Mj- 
For directed networks, recent approaches are proposed to 
ensure this convergence [19]. 

The random walk gives rise to a geometric diffusion 
with an associated distance between nodes i and j de- 
fined as [17|]: 4(t) = Y.k>o ^^^i^^^^^, where the term 
/i^ is supposed to compensate for discrepancies in local 
densities. By construction, distance between nodes is 
strongly ruled by the connectivity of the graph, and it 
takes small values if nodes are connected by many paths. 
Considering the spectral representations of matrix P, one 
has a set of eigenvalues |Ao| > |Ai| > ... > |AAr-i| 
and eigenvectors ipk and V^/e such that (^^P = Xk^J. 
and Pil)k = Xk'^k- The diffusion distance can be writ- 
ten as: dj.{t) = T.k>i^fi^k{i) - ^k{j)f where V^/e(j) 
denotes the component j of eigenvector k. The diffu- 
sion distance can be approximated (note that ipQ = ji" 
and t/^o is a constant vector) to a relative precision us- 
ing the first /3 nontrivial eigenvectors and eigenvalues: 
4(^) - ELi^n (V^n(0 - V^n(j))'. This is equivalent 
to embed the graph in a low dimensional space R^, con- 
verting the diffusion distance between nodes of the graph 
into Euclidean distance in [I^. This approach has 
the main advantage that it defines a meaningful represen- 
tation of the graph and it leads to explicitly define a dis- 
tance metric on the space that reflects the connectiv- 
ity of the network. Graph modules are then extracted by 
a /c-means clustering algorithm in the embedding space. 
The algorithm starts with a random assignment of Nm 



cluster center. Then, the partition is updated by repeat- 
ing the following steps: i) each point is assigned to the 
nearest center; ii) the new geometric center of each clus- 
ter are recomputed. The algorithm stops when the as- 
signments of nodes are no longer changed. The algorithm 
reaches convergence in very few iterations after the initial 
generations, but it only guarantees a convergence towards 
a local optimum. In fact, the /c-means is sensitive to the 
initial choice of centroids. To overcome this drawback, we 
run the algorithm several times (500) for different num- 
ber of modules {Nm = 2, . . . , 15) and return the partition 
yielding the largest modularity Q. The modularity Q{S) 
for a given partition 5 of a weighted network is defined 

as [^: Q{S) = ^^^-^ - (||;)^j , where L is the total 
weight of all connections in the network, Ig is the weight 
of links between vertices in module 5, and dg is the total 
weight of links in module s. To select the optimal number 
of diffusion coordinates, we have repeated the clustering 
algorithm for all possible values of /3, and looked for the 
largest value of Q. 
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FIG. 2. (a) Brain sites belonging to each functional brain 
module were colored and superimposed onto an anatomical 
image (there is no color correspondence between the modules 
of patients and those of controls). For the sake of clarity, iso- 
lated nodes were colored in black. Topographic distribution 
and difference maps (p- values) of participation coefficient (b) 
and within-module degree (c) for real and equivalent random 
graphs. 

The method has been tested on synthetic networks 
[20], and used to find modules in brain networks. The 
presence of modules is actually confirmed by the high val- 
ues of Q obtained for the optimal partition of the brain 
networks of both control subjects and patients. The main 
result is that brain networks of epileptic patients have 
a structure consisting in a smaller number of modules 
{Nm = 8.30) with respect to that of control subjects 
{Nm = 13.7). To assess the stability of the partition 
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structure across healthy subjects and patients we used 



the adjusted Rand index J [2l|, which yields a normal- 



ized value between (if the two partitions are randomly 
drawn) and 1 (for an identical partition structure). The 
values of J indicate a high stability of the partition struc- 
ture across all the epileptic patients (J = 0.733) and cer- 
tain variability in the modular structure of control sub- 
jects (J = 0.479). 

Fig.[2fa) illustrates the spatial distribution of the mod- 
ules obtained for healthy subjects and patients. De- 
spite the differences, some modules fit well some known 
brain regions including occipital, parietal, and frontal ar- 
eas. Although a one-to-one assignment of anatomical 
brain areas to the retrieved modules is difficult to de- 
fine, modules assignment provides the basis for the anal- 
ysis of nodes according to their patterns of intra- and 
inter- modules connections [5]. The within- module de- 
gree z-score measures how well connected the node i is 
to other nodes in the same module, and is defined as 
-; where kis is total weight of links of node 



i to other nodes in its module 5, while kg and (Jk^ are 
average node strength and standard deviation of nodes 
in s. Thus node i will display a large value of Zi if it has 
a large number of intra-modular connections relative to 
other nodes in the same module. The participation coef- 
ficient pci = 1 — Yl^^i quantifies instead to which 
extent a node i is connected to different modules. This 
coefficient takes values of zero if a node has most of its 
links exclusively with other nodes of its module, and 1 if 
they are distributed among different modules. 

In Figs.[2fb-c) we report respectively the spatial distri- 
bution of the node participation coefficient pc and of the 
z-score for control subjects and patients. The difference 
map for the z-score reveals very small differences in the 
way the nodes are connected to other nodes in the same 
module. Conversely, the distribution of the participa- 
tion coefficient pc in control subjects strongly differ from 
that obtained for epileptic patients where nodes (spe- 
cially those of the right centro-parietal areas) participate 
with their links in several modules. The same modular 
partition did not reveal significant differences between 
the equivalent random configurations of both groups [22]. 
This is a remarkable result as it supports the hypothesis 
that normal and pathological brain dynamics are charac- 
terized by different functional connectivity patterns [23] . 

In conclusion, in this paper we have addressed a fun- 
damental problem in neuroscience: characterizing the 
connectivity structure of functional networks associated 
to normal and pathological neural dynamics. From the 
analysis of spontaneous brain activities at rest, we found 
that the architecture of functional networks extracted 
from epileptic patients differs from that of healthy sub- 
jects. Interestingly, we identified a non-random modular 
structure of brain networks. Modularity analysis revealed 
that nodes of epileptic brain networks abnormally link 



different functional modules in the network. The con- 
nectivity of brain activities at the extended a-band of 
epileptic patients might play a putative role in the emer- 
gence of absence seizures. This leads up to the need of 
more refined studies, as in Ref [24j, to address the role of 
this architecture in the absences seizures. 

A modular description of brain networks might pro- 
vide, more in general, meaningful insights into the func- 
tional organization of brain activities recorded with oth- 
ers neuroimaging techniques (EEG, MEG or fMRI) dur- 
ing diverse cognitive or pathological states. Applied to 
other multivariate data (e.g., financial or epidemiologi- 
cal time series), our approach could provide new insights 
into the network structure of spatially extended systems. 
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